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ABSTRACT 

We have constructed an analytical model of AGN feedback and studied its implications 
for elliptical galaxies and galaxy clusters. The results show that momentum injection 
above a critical value will eject material from low mass elliptical galaxies, and leads to 
an X-ray luminosity, Xx, that is oc a^~^^ , depending on the AGN fuelling mechanism, 
where a is the velocity dispersion of the hot gas. This result agrees well with both 
observations and semi-analytic models. In more massive ellipticals and clusters, AGN 
outflows quickly become buoyancy-dominated. This necessarily means that heating 
by a central cluster AGN redistributes the intracluster medium (ICM) such that the 
mass of hot gas, within the cooling radius, should be oc Lx(< 7'cooi)/[.9(?'cooi)o'], where 
^(''cooi) is the gravitational acceleration at the cooling radius. This prediction is con- 
firmed using observations of seven clusters. The same mechanism also defines a critical 
ICM cooling time of ~ 0.5 Gyr, which is in reasonable agreement with recent obser- 
vations showing that star formation and AGN activity are triggered below a universal 
cooling time threshold. 
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1 INTRODUCTION 



It is well known that e lliptical galaxies are commonly the 
hosts of powerful AGN ()McLure et al.ll2004h . These sources 
give rise to lobes of radio emission embedded in the X-ray 



ating clusters of galaxies (e.g. iBurns 1990; Best et all 


2005 


Birzan et al.l 2004: Dunn et al.l 


2005,: Raffertv et al. 


200e 


BestI 20071: Nulsen et al. 2007: Diehl & Statlei 20081). 



Recent observational studies of AGN in elliptical galax- 
ies, and galaxy clusters, also suggest that AGN activity 
is related to the thermal state of its environment. Sys- 
tems with short radiative cooling times, or a low central 
entropy, have been shown to be more likely to host ac- 
tive star formation, optical linc-cmission and jet-produ cing 
AGN fe.g. [Burns 1990: Craw ford ct al. 1999; Raffcrty e t"al] 
l2008l : ICavagnolo et all l2008l : iMittal et al.li2008t ). This sug- 
gests that AGN activity is part of a feedback loop that is 
likely to have important consequences for its environment. 

The theoretical evidence for the impor tance of AGN 
feedba ck h as also continued to m ount. Ta bor fc BinnevI 
l|l993h and iBinnev fc Taboj l| 19951 ) employed AGN feed- 
back to prevent catastrophic radiative cooling of the X- 
ray emitting gas that surrounds elliptical galaxies. If un- 
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opposed, radiative cooling would lead to the deposition of 
large quantities of cold gas in the central galaxy, and thus 
cause the galaxy to grow rapidly. Since this appears not 
to be the case, it was assumed that there must be a heat 
source that balances the cooling. In a different approach to 
a similar problem, s emi-analytic inodels of galaxy formation 
and evolution (e.g . iBenson et al.l l2003l : ICroton et al.l l2006l : 
iBower et al.l l200d ) also demonstrated that AGN feedback 
was necessary to reproduce the observed galaxy luminos- 
ity function at large masses, and balance the large radiative 
cooling rates in clusters of galaxies. In the absence of AGN 
feedback, the models produced too many massive galaxies. 

Semi-analytic models (e.g. Benson et al.l 20031), hy- 
drodynamic simulations (e.g. 
Gael j [l989l: IPope et all l2005l) 



Bregman fc DavidI iTgisI : 
and observations (e.j 



Voigt fc FabianI |2004| : IPope et al.) |2006| ) have all shown 



that thermal conduction cannot provide the required en- 
ergy transfer rates to balance cooling in the majority 
of cases. However, numerous hydrodynamical simulations 
have shown that a combination of AGN heating and ther- 
mal condu ction can maintain the system in a steady- 
state (e.g. Ruszkowski fc BegelmarJ 20o3: Briiggeni l2003l : 



Guo et all ( 2008h for a stability analysis. 



igge 

Hoeft fc Briiggenll2004l : lRovchowdhurv et al.ll2004l ). see also 



Recent theoretical efforts have focussed on trying to un- 
derstand the X-ray luminosity/temperature (Lx ~ T) rela- 
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tion in elliptical galaxies an d galaxy clusters. The results 
of sem i-analytic models (e.g. iBower et al.l 20081: Dave et al] 
I2OO8I ) and hydrodynamic models (e.g. lPuchwein et al.ll2008l l 
have shown that AGN feedback does indeed play a signifi- 
cant role in governing the scaling of this relation. 

The aim of this article is to describe the impact on el- 
liptical galaxies of momentum and thermal energy injection 
by a central AGN. There are critical thresholds above which 
each energy component will eject material from the galaxy, 
and so each mechanism leads to clear observable signatures 
which can be used to distinguish between them. This study 
is restricted to the study of galaxies where the X-ray emis- 
sion is dominated by the gaseous atmosp here, rather than 
stellar sources (e.g. lO'SuUivan eraLll200lh . 

To this end, we have employed an analytical model of 
AGN feedback in order to better understand, and compli- 
ment, the results that have emerged from semi-analytic and 
hydrodynamic simulations, as well as observations. 

The article is arranged as follows: in section 2 we in- 
troduce the basic model describing the efi'ect of energy and 
momentum injection into a hot gaseous atmosphere. Section 
3 contains the results for the limiting cases of momentum 
dominated and thermally dominated outflows and discusses 
the implications for the Lx — o" relation in elliptical galaxies. 
Note, we have chosen to study the Lx — cr relation rather 
than Z/x — T, because the gravitational potential can be 
written in terms of a which we assume to be unaltered by 
feedback. In section 4, the theory is applied directly to obser- 
vations of galaxy clusters. The main findings are summarised 
in section 5. 



2 THE MODEL 

The aim of this section is to develop a simple, but use- 
ful, analytical formulation that can describe the general im- 
pact of AGN feedback on the Lx — cr relation of elliptical 
galaxies. For this purpose the elliptical g alaxy is approxi- 
mated as a singular isothermal sphere fe.g. lKinal2003l . [20051 : 
[Murray et al. 2005). This is because the setup has proved 
to be a successful starting point in describing phenomena 
such as the black hole mass/stellar velocity dispersion rela- 
tion. The results are used to illustrate general behavioural 
trends, and are not intended to be a complete description 
of all the physical processes involved. However, in Section 
4, which focuses on galaxy clusters, observational data are 
used to calculate more accurate gravitational potentials and 
gas masses. This allows specific predictions of the model to 
be tested rigorously on individual systems. 

For the singular isothermal sphere, the mass of gas 
within a radius R, is 



(1) 



where /g « 0.17 (|Spergel et al.ll2007l ) is the gas fraction, a 
is the velocity dispersion, R is the radius, and G is New- 
ton's gravitational constant. In this model, the AGN resides 
at the centre of the potential, and produces outflows that 
propagate radially outwards. The mass swept up by out- 
flows subtending a solid opening angle, Q, is therefore 



so that M{R) = r2MT(-R)/(47r). As long as the outflow is 
active, new material from the gaseous atmosphere will be 
prevented from entering this solid angle through the sides. 
Therefore, equation© gives a reasonable estimate of the 
mass collected by the outflow. The darkmatter gravitational 
potential is assumed to be unaffected by the movement of 
the gas. 

Considering momentum injection at a rate of P, New- 
ton's 2nd law gives 



d[M{R)R] 
dt 



= P + g{R)M{R), 



(3) 



where g{R) is the gravitational acceleration. If P ^ 
~g{R)M(R), momentum injection will overcome the gravity 
of the galaxy and the AGN will induce a bulk outflow of hot 
material collected from the galaxy's atmosphere. 

However, the outfiow must do work against both gravity 
and against the pressure of the gaseous atmosphere it is 
lifting. This means that a term, {Q,/A-k)A-rR^p^^\,, must be 
subtracted from the right hand side of equation ((3)l, where 
Pamb is the ambient pressure. Using the ideal gas equation, 
the pressure term can be written 



^ )47rPVmb(P) 



f\,.R^ 
47r I 



p{R)kbT _ / Q.\ 2/gcr'* 
(4) 

where ky^T/finip = a^. 

The AGN also injects thermal energy, at a rate -Ethcrmai, 
that will exert an outward pressure on the ambient gas. As- 
suming that the infiation is reversible 



^thermal = ( ^ R' Rp,^^{R) , (5) 

7 — 1 \ 47r / 

where 7 is the adiabatic index of the injected gas. Therefore, 
we must add a term 

^ I /I d2 I r,\ 7 ~ 1 -Ethcrmai , „s 

— HttP Pamb(P) = ;i; , (6) 

to the right hand side of equation (IS}. (This ter m should be 
comp ared with the LEddington/c term used in iKind l2003l . 
l2005l . to model the outflow momentum.) The completed 
equation of motion can be written 

MM . P+,(P)M(P)- (^) ^+1^^. 

(7) 



3 RESULTS 

In this section we investigate the behaviour of equation 
(O in two limits: momentum-dominated AGN heating, and 
thermally-dominated AGN heating. 



3.1 Kinetic feedback 

In the limit that the outflow is completely momentum 
driven, the critical value of the momentum injection rate 
for the singular isothermal sphere is 



R 



(2) 



P„. = -ff(P)A/(P) + ( j^) ^ 



6^ 
G 



(8) 
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In general, it is more intuitive to consider the energy injec- 
tion rate rather than the momentum flux. Therefore, equa- 
tion ([8]) shall be re-written accordingly. The momentum in- 
jection rate is related to the kinetic energy injection rate 
by 



2mou 



(9) 



where rhout is the mass flow rate of the jet. 

Using equation ^ with equation ((HJ and rearranging 
for the critical luminosity gives 



nctic.crit — 



which, in scaled units, is 



^kinetic, crit 



2 X 10' 



[ 0.02 



(10) 



lOMoyr- 



200 kmh 



erg s 



(11) 



The reason Q/Att = 0.02 has been chosen will become ap- 
parent in section 4, where the value is derived from observa- 
tions. At present, it is sufficient to note that f2 is an effective 
opening angle, rather than the actual opening angle of the 
outflow. It is the solid angle within which the outflow in- 
teracts with ambient material, and is therefore larger than 
the typical value derv ied from radio observations (e.g. see 
iBinnev fc Taboilll995l . and references therein). 

The luminosities hinted at in equation (|ll|l mean that 
it is plausible to expect kinetic AGN feedback to eject hot 
material from the atmosphere of the galaxy. This will have 
two effects: it will reduce the X-ray luminosity of the galaxy, 
Lx, and reduce the gas fraction. 

In clusters (with large a), the critical jet power will 
exceed the cluster X-ray luminosity, and would require im- 
plausibly large jet powers to eject mass. This means that the 
AGN jets are unlikely to significantly affect the Lx oc a^'* 
relation (jMahdavi fc GcUcr 2001) in clusters. Of course, this 
does not mean that they do not play a significant role in 
heating the inner regions of such objects. 

To make a realistic estimate of the Lx — cr relation ex- 
pected as a consequence of momentum -dominated feedback, 
we follow the method o f iBower et aD (,2008). Their prescrip- 
tion for the effect of the AGN activity was to remove mass 
from the galaxy at a rate governed by the difference be- 
tween the heating and cooling rates: m — 2{E — Lx)/cr^. 
Our model is expected to behave in a similar way, but with 
the reference point being the critical jet power defined in 
equation (|ll|l 



2(-Ekin 



(12) 

Equation (|12p tells us that the mass-loss from the galaxy will 
effectively stop when the average feedback heating rate be- 
comes comparable with the critical jet power for that grav- 
itational potential. 

The complete derivation requries a functional form for 
L^kinetic- [Sower et af. (2008) used an energy injection rate 
that is proportional to the mass cooling rate, Afcooi, of the 
X-ray emitting atmosphere around the galaxy. We will follow 
their example as an illustrative case. Therefore 



-Ekinotic ~ r]SMcoolC 



(13) 



where the factor S accounts for the fraction of cooling ma- 
terial that reaches the black hole. Here, S is also used to 
account for the reduced cooling rate compared to the classi- 
cal value. Given that cooling in such systems appears to be 
reduced by a factor of 10, and assuming that only one tenth 
of the infiowing material reaches the black hole, S ~ 0.01. r; 
is the accretion efficiency, and is the speed of light. The 
classical mass cooling rate can be written in terms of the 
X-ray luminosity of the hot g 

M - 

JWcool — TT, 

5 

so that the kinetic energy injection rate is 

„2 



(14) 



-Bkin 



: OAr]SLx- 



(15) 



Substituting equation (|15|l into equation (I12p shows that the 
mass-loss rate from the galaxy will tend to zero as 



2.5- 



^kinetic, crit 10 
OC a . 



(16) 



In scaled units, this means we expect the X-ray luminosity 
of an elliptical galaxy to be approximately 



Lx ~ 2 X 10' 



V 0.02 



_5_ 

ooT 



friout 



200 km s 



lOMoyr-i 

10 

ergs 



(17) 



where an accretion efficiency of r; — 0.1 has been 
assumed. These values agree reasonably well with the 
observed X-ray lum inosities of elliptical galaxies (e.g. 
iFukazawa et al.ll2006l '). The resuh agrees well wi th the ob- 
servations presented bv lMahdavi fc Gellerl |20o3) who also 



found L 



X oc g 



and is consistent with the relation found 



bv lO'SuUivan et al.1 (|2003|). The normalisation of the rela- 
tion found bv lMahdavi fc Gelled | 200 ll ) seems to be slightly 
higher than equation (|17p . possibly suggesting that 5 ~ 10~^ 
would be a better description of the physics governing 
the accretion rate. This would guarantee a low Eddington 
fraction accretion rate which would produce outflows (e.g. 
iMaccarone et al1l2003l ) and not quasars. 

The scaling relation in equation (|17|l is steeper t han the 
results of numerical simulations performed by Dave et al.l 
(|2008l ) who found Lx oc a^'^. However, our model would 
also produce Lx oc a* if the AGN fuelling rate was propor- 
tional to Lx rather than Lx/cr^. In both observations and 
simulations, there is signiflcant scatter in these relationship. 
Presumeably this could be caused by the different histo- 
ries of individual galaxies (mergers and ages) and also the 
stochastic nature of AGN feedback. 



3.2 Thermal energy injection 

It is also important to consider the role of thermal energy 
injection. According to equation (O, the effect of thermal 
injection should be quite different to pure momentum injec- 
tion. This leads to several predictions which can be tested 
against existing observations. 

A variant of this mechanism is generally considered to 
lead to the black hole mass - bul ge stellar velocity disper- 
sion relation (e.g. IPuchwein et al.1 ,2008). However, in their 
approach the radiative power output of the AGN accretion 
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at a high Eddington ratio couples to the ambient gas. We will 
not be considering effects of radiation output from the AGN 
since it is assumed here that the AGN simply injects ther- 
mal energy that does work against its surroundings, rather 
than assuming anything more complex. 

Using equation ((TJ, the critical energy injection rate 
for the thermal gas to eject material from the gravitational 
potential is of the order 



thermal, crit 



4-K G 



where we have set _R ~ ct. In scaled units this gives 



thermal, crit 



10^ 



V 0.02 



200 km s- 



erg s 



(18) 



(19) 



In reality, the AGN injects both momentum and thermal 
energy (as well as the radiation that we are not considering) , 
rather than one or the other. The behaviour of the system 
will therefore depend on the relative proportions. 

For simplicity, the thermal energy injection rate can be 
written as a fraction, x, of the total energy injection rate; 
^^thermai = X^totai . Therefore, the critical threshold for ther- 
mal energy injection, in terms of the total energy injection 
rate, is Stotai.crit = £^thermai,crit/x- It is probably safe to as- 
sume that initially, at least, outflows are momentum domi- 
nated, since they are highly anisotropic. This suggests that 
X ^ 1. In the next section we will show that, above a critical 
jet length, AGN outflows will become buoyancy dominated 
so that X ~* 1 for large radii. However, for the remainder 
of this section it is assumed that the relative fractions of 
kinetic and thermal energy in the outflow remain constant, 
and that the outflows are momentum dominated. 

Combining the critical luminosities for both thermal 
and kinetic limits we flnd that the critical energy injection 
rate can be expressed as a broken power-law that scales as 
cr* at lower masses, and at higher masses. The location of 
the break depends on the value of x; (crbroak/200 km s~^) ~ 
(2x)~^''"^. So, as an example, if x = 0.1, the break will oc- 
cur at CT ~ 340kms~^. Above this value of a, the criterion 
for mass ejection from the galaxy follows the thermal limit, 
while below this the mass ejection criterion is governed by 
momentum injection. Another way of expressing this is that 
thermal energy injection has much more important conse- 
quences for the galaxy properties at high masses, while mo- 
mentum injection is more important at lower masses. This 
result is shown in figure 1. 

The critical energy injection rate at cr ~ 340kms~^ is 
^ 10^** ergs~^ which is well below the Eddington limit for a 
10* Mq black hole. This strongly suggests that thermal feed- 
back from black holes could have an effect on the Lx — cr rela- 
tion in massive elliptical galaxies and groups. However, for a 
cluster gravitational potential, a > 500kms~^. The critical 
energy injection rate a,t a — 500kms~^, is ~ lO^'^ergs"^, 
which is comparable to the Eddington limit for a 10^ Mq 
black hole. Furthermore, depending on the value of 5 the 
actual power available from accretion might be much lower 
than this. Consequently, there should be a second break 
in the Lx — a relation at cr ~ 500kms~^. Assuming that 
— k]jT/fj,mp, a — 500kms~^ corresponds to a tempera- 
ture of ~ 1 — 2 keV, whi ch is approximately where a break is 
observed (e.g. IXue fc Wu 2000'). This, therefore, marks the 
transition above which black holes are unlikely to affect the 



10" - 



10" - 



100 



velocity dispersion (km s ^) 



Figure 1. The critical energy injection rate for thermal (dashed 
line) and momentum (dotted line) injection. In this example it 
was assumed that the thermal fraction, x = 0.1 so that the critical 
values cross at roughly 340kms~^. At low a, the lower critical 
energy injection limit is provided by the momentum criterion, 
while at high a the lower critical injection limit is provided by 
thermal energy criterion. 



global Lx — cr relation. Nevertheless, black holes may still 
be effective in heating the central regions of clusters. 

If the accretion rate onto the black hole is oc Ly^/a^ , as 
suggested in the previous example, thermal AGN feedback 
should lead to Lx oc in higher mass ellipticals and groups. 
Alternatively, if the black hole accretion rate is oc Lx, the 
model would lead to Lx oc . 

So, in summary, the model predicts Lx oc cr*"'^'^ for 
lower masses elliptical galaxies, Lx oc a^~'^ for higher mass 
ellipticals and groups, depending on the AGN fuelling mech- 
anism. However, given that the relation for c lusters can 
be written Lx oc a*'*' l|Mahdavi fc Gelled [iooH ) and there 
is a noticeable b reak in the scaling at ~ 1 — 2 keV (e.g. 
IXue fc Wul I2OOOI ') (or a ~ 500 km s"^ ), it seems that an 
AGN fuelling rate that is oc Lx/f^ can best explain the 
observed Lx — f relation across the mass range considered 
here. 



4 APPLICATION TO GALAXY CLUSTERS 

Having considered the implications of momentum and ther- 
mal energy injection on elliptical galaxies, the model is ap- 
plied to clusters of galaxies. Because of the stronger gravity, 
these systems behave rather differently to elliptical galax- 
ies. This leads to several observable consequences which are 
discussed below. In addition, the benefit of using the singu- 
lar isothermal sphere for elliptical galaxies was to construct 
simple scaling laws across the population of galaxies that re- 
sult from different properties of AGN feedback. However, the 
previous section demonstrated that AGN feedback is unable 
to alter the global Lx — o" relation in clusters, so using the 
same general approach would not yield interesting results 
for clusters. Instead, it seems more appropriate to focus on 
the central regions of clusters, where the AGN is likely to 
have an impact and so attempt to discern any observable 



1000 
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signatures of AGN feedback. To do this accurately requires 
observations of the temperature and density profiles of the 
central ICM in a sample of clusters. This is done in Section 
4.2. Firstly, though, there is brief discussion of the general 
effect of the cluster potential on AGN outflows. 



4.1 Buoyancy, work and bubbles in the 
intracluster medium 

Momentum and thermal energy fractions govern various 
properties of the outflow, including its morphology. These 
fractions are not constant, but vary as the outflow propa- 
gates through its surroundings. For example, thermal energy 
injection produces buoyant material, and the buoyancy force 
is also a function of the local gravitational acceleration. This 
suggests that identical jets will appear morphologically dif- 
ferent in different environments. In particular, buoyancy is 
more important in more massive ellipticals, and in clusters, 
due to the stronger gravity in such enviroments. It, therefore, 
stands to reason that we might expect outflows to appear 
more bubble-like (buoyant), and less elongated, in clusters. 
This possibility is investigated below. 

For a given thermal energy injection rate, it is pos- 
sible to calculate a buoyancy flux. An outflow of veloc- 
ity, w, and width, D has a buoyancy flux defined as (e.g. 
ICarlotti fc Hunj|2005l . and references therein) 



„ _ Ap nD^ 4 -3 

B = gw cm s , 

P 4 



(20) 



where g is the gravitational acceleration, p is the ambient 
density, and Ap is the density difference between the outflow 
and the ambient material. The buoyancy flux can also be 
written in t erms the heating rate that caused the density 
change (e.g. ICarlotti fc HuntJbOOSh 

7-1 ggthcn^al ^^4^-3^ (21) 
7 Pamb 

where Pamb is the pressure of the a mbient gas. Similarly, the 
momentum flux is deflned as (e.g. ICarlotti fc^ Huntll2005l ) 



U 



^2^')^2^ = :^cmS-^ (22) 
p J 4 p ^ ^ 

where poutflow/p is the ratio of the density in the jet to the 
ambient density. 

It is useful to construct a characteristic length-scale 
from the momentum and buoyancy fluxes. This length-scale 
marks the transition from a momentum-dominated flow to 
a buoyancy-dominated flow and is deflned as 

[/3/4 



(23) 



Using the fact that P = (2mouti5kinctic)^^'^, the length-scale 
can be written 



/ ~ 24 kpc 



-3/4 



JTlout 



3/8 



lOMgyr-i J \W^'^ ergs-i 



3/8 



10- 



' g cm" 



Pamb 



1/2 / p ^ -1/2 

/ -^thermal 



10^^" ergcm~3y \10"'^ergs~i 



9 



-1/2 



10~® cms~2 

We have assumed a uniform density for the environment. 



and constant gravity, but nevertheless, the critical length 
above which the jet dynamics become dominated by buoy- 
ancy is likely to be ~ 20 kpc. 

The parameters that control this length-scale show that 
an outflow will become buoyancy-dominated closer to its 
origin in a denser environment with stronger gravity. This 
confirms the earlier statement that, for a given energy in- 
jection rate, bubble-like outflows will be more prevalent in 
the centres of clusters of galaxies than in low mass elliptical 
galaxies. It also means that more of the injected energy is 
trapped within the central regions, providing more efficient 
heating in a cluster environment. 

In clusters, it is not possible to relate many of the prop- 
erties, such as density and pressure, to a. As a result, it is 
not particularly useful to represent equation (|24p in terms of 
the velocity dispersion. In addition, there is no single value 
of a when the jets become buoyant, because this will de- 
pend on the kinetic and thermal energy injection rates and 
the other cluster properties. 

Another possible explanation for the presence of AGN- 
blown bubble-like structures in clusters is that jets never 
exceed the critical energy injection to eject mass from a 
cluster. Therefore, the jet will always be constrained by the 
ambient medium. In contrast, it is possible for the AGN 
to exceed the critical luminosity in an elliptical galaxy, in 
which case the jet is no longer constrained. It is possible 
that this effect plays so me role in the apparen t FRI/FRII di- 
chotomy. For example, iKaiser fc Besd |2007i) suggested that 
FRIs were those objects which are disrupted within the core 
region the surrounding atmosphere, while FRIIs were jets 
which were sufficiently powerful to penetrate into the re- 
gion outside the core where the ambient pressure drops off 
rapidly. Note, however, other effects such as black ho le spin 
are a lso a possible explanation for the dichotomy (e.g. lMeieil 
Il999h . 

It is also necessary to show that the term i?thcrmai/A 
term in equation (O is equivalent to the work done by 
bouyancy against the surrounding gas. The rate of change 
of the buoyancy force, according to Archimedes' principle, 
can be written 



F = 



■iivpgR^R, 



(25) 



where the density of the outflow is assumed to be much less 
than the material it has displaced. Using equation (|2ip the 
rate of change buoyancy force in terms of the heating rate, 

-S'thcrmali is 



F = 



7 - 1 pmpffi^thcrmal 

7 kbT 



(26) 



Equations (|25p and (|26p can then be rearranged to show 
that the two forms are equivalent 



7 



R 



^ ]47rpamb-R^ 



(27) 



This means that buoyancy has the same effect as thermal 
energy injection doing work against its surroundings. Fur- 
thermore, since the outflows become buoyancy-dominated 
on length-scales greater than I, this means that the initial 
(2ifc)omentum flux is converted into thermal energy. The total 
power must be conserved, so for R > I, the outflow can be 
adequately described by ^thermal ~ Stotai, where Btotai is 
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the (constant) total energy injection rate of the outflow. As 
a result, the factor x, introduced earlier, is a strong function 
of radius and becomes close to unity outside R = I. 

Equation (|24p suggests that I is likely to be significantly 
less than the cooling radius in a typical cluster (~ 100 kpc). 
Therefore, we argue that on spatial scales comparable with 
the cooling radius, it is sufficient to model heating by an 
AGN at the cluster centre, purely in terms of the injection of 
thermal energy within a given solid angle. To understand the 
dynamical structures on smaller scales we must consider the 
jet momentum in detail, but on intermediate spatial scales 
it seems that the energy injection can be treated as purely 
thermal. We will next consider the physical implications of 
thermal heating in galaxy clusters. 

4.2 Observational consequences of buoyancy in 
cluster cores 

In clusters, the AGN cannot eject gas from the cluster due 
to the sheer weight of material residing in the gravitational 
potential. Ifowever, there is still a region near the clus- 
ter centre where the AGN can affect the distribution of 
the X-ray emitting material. Since we are considering clus- 
ters, buoyancy is important and the energy injection can 
be thought of as purely thermal. Therefore, the appropriate 
critical energy injection rate is the thermal injection limit; 
Ecrit ~ -{15/4)(TM{R)g{R). As usual, M{R) is the mass of 
the material swept up within R and within the solid angle 
of the outflow; M{R) = Qj^r^p{r)dr = MfT(7?)/(47r). a 
is the velocity dispersion of the gravitational potential, and 
g{R) is the gravitational acceleration at R. In general, clus- 
ters are close to hydrostatic equilibrium so that the gravita- 
tional acceleration can be calculated from the pressure gra- 
dient and the gas density; g{R) « Vp{R)/p{R). The multi- 
plicative factor 15/4 comes from (3/2) x 7/(7 — 1), assuming 
that 7 = 5/3. Finally, the factor 3/2 is from the contribution 
of the ambient pressure to the total work required for the 
thermal outflow to eject material, see equation (0. 

For energy injection rates less than this value, material 
can build up in the central regions of the cluster, while for 
larger energy injection rates, mass is ejected from the central 
regions. This is essentially the same as the explanation for 
mass loss from the elliptical galaxies, except that in clusters 
mass is ejected from the central regions, but still retained by 
the cluster as a whole. Of course, this implicitly assumes that 
some of the energy injected by the AGN is distributed over 
47r steradians. This could occur through shocks generated 
in th e early stages of outflow and subsequently sound waves 



Voit fc Donahue! l 2005h . or bulk motions in the ICM 
Heinz et al. I bood r though these mechanisms cannot 



(e-g- 
(e-g- 

be accounted for in the current model. 

The general argument outlined above suggests that 
there will be a radius, governed by the energy injection rate, 
within which the weight of X-ray emitting material is gov- 
erned by AGN feedback. This is important, since by con- 
trolling the mass of gas in the central regions of the cluster 
the AGN also regulates the cooling rate, and in doing so its 
fuel own supply, and so on. 

For clusters, the obvious region of importance is en- 
closed by the cooling radius, since it contains the material 
that could cool within the age of the Universe. Therefore, 
if AGN feedback balances the radiative losses within the 



cooling radius, E — Lx(< r'cooi), then, according to the 
model, the heating should also redistribute the material in 
the centre of the cluster. This interplay should have reached 
a quasi-equilibrium when 



E = Lx{< rcooi) = 

(t) 



CrA/(rcool)ff(rcool) ~ 



-47r 



Vp(rcooi) 



rp{r)dr. (28) 



P(rcool) 

In other words, the total mass of X-ray emitting material 
within Tcooi, multiplied by the gravitational acceleration at 
''cool should be 

15/ V 



A'/T(rcool)g(''cool) 



(29) 



To be strictly correct, the AGN removes a fraction, S7/47r, 
of this material, so the remaining weight must be —(1 — 
r2/47r)MT(rcooi)5(7-cooi). 

For comparison, AGN energy injection rates can be es- 
timated from observations o f bubble enthalpies divided by a 
characteristic timescale (e .g. iBirzan et al ] |2004| :Id unn et al.l 
I2OO5I : iRaffertv et aLll2006l). Although these values are some- 
what uncertain fsee lNusser et al]l2006l : [Binnev et al ] |2007l . 
for example), it seems that more powerful AGN outbursts 
do seem to be observed in more luminous clusters as well 
as those systems with shorter central cooling times. Thus, 
AGN outbursts in clusters are capable of injecting energy at 
the rate required by the model. 

Equation (|29p has been veri fied by examining the sam- 
ple of seven clusters discussed in lPope et al.l (|2006l 'l and can 
be seen in figure 2. The sample provides a readily available 
data set with all the physical parameters required to inves- 
tigate the model predictions. The main reason for this is 
that the clusters are nearby so that the temperature and 
the density profiles are fairly well determined to within a 
few kiloparsecs of the cluster centre. Accurate temperature 
and density profiles permit reasonable estimates of the total 
gas mass, and gravitational acceleration, thus increasing the 
possibility of detecting subtle variations due to AGN feed- 
back. Each member of the sample is known to host a central 
AGN, thus making an ideal testbed for the model. Due to 
the paucity of very high X-ray luminosity clusters, future 
work will focus on extending the sample predominantly at 
low X-ray luminosities, to investigate and refine the model. 

The best-fit in figure 2 implies a typical solid opening 
angle of n/47r ~ 0.028, which informed the value of the 
opening angle used throughout the article. Noting that this 
value is twice the effective solid opening angle of the outflow, 
the effective half opening angle of a single outflow must be 
~ 14 °. This is a factor of ~ 3 greater than the typ ical values 
found by radio astronomers (e.g. see lBinn"ev" fc Tabor 199^, 
and references therein). Possible explanations for this are 
discussed below. 

In the simplest case, a buoyant bubble inflated by an 
AGN would expand adiabatically to maintain pressure equi- 
librium with its surroundings. According to this model, the 
bubble radius sh ould behave in the following manner (e.g. 
iDiehl et al.ll2008h 



Pa,mh (R) 



PO 



i/3r 



(30) 



where r\^{R) is the bubble radius as a function of radial dis- 
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placement from its original location. rb,o is the initial radius 
of the bubble, po is the pressure at this location, and Pamb(^) 
is the ambient pressure as a function of radius. F is the poly- 
tropic index of material inside the bubble. The half opening 
angle, 9, can then be calculated from tan{9) — drb(ii)/d_R. 
Substituting suitable values for rb,o, and F = 4/3 and us- 
ing a generic cluster pressure profile suggests that the adia- 
batic expansion of bubbles would not result in large enough 
opening angles to explain the b est-fit in fig ure 2. In a more 
realist ic model. [Pavlovski et al.l (2007) and IPavlovski et al.l 
l|2008t) showed that the Kutta-Zhukovsky force causes bub- 
bles to expand further than would be expected from pressure 
equilibrium arguments. 

Probably the most plausible explanation is that the ef- 
fective opening angle derived from figure 2, refers to the 
entrainment of material from the boundary layer between 
the outflow and the ambient medium. More specifically, en- 
training material along the length of the outflow, as well as 
through the head of the outflow. Consequently, the outflow 
can collect material from within a larger s olid angle than 
is carved out by th e main flow itself, (see iBicknelll 1 1994 
iLaing fc Bridi3l2002l . for example). 

Generally, 9 is re lated to the entrainment coefficient, 
a, by Q = tan(^) fe.g. lTurneijll986l ). In a uniform ambient 
medium, and in the Boussinesq limit, the typical entrain- 
ment coefficient for a momentum driven outflo w is q ~ 0.05 , 
but a ~ 0.08 for a buoyancy driven flow (see lTurneij|l98(j . 
and references therein) . The exact values depend slightly on 
whether the cross-section of the outflow is modelled as a 
Gaussian or top-hat function, but the general result is that 
entrainment coefiicients are larger for buoyancy driven flows. 

The entrainment coefficient also varies with the density 
contrast between the outflow and the ambient m aterial, and 
is a f unction of the Mach number of the flow. iDe Youn3 
(| 19931 ) suggests the opening angle can be written as 



tan(6l) ~ 0.08 



Poutfiow 



1/2 



(31) 



where p is the ambient density, poutflow is the density of the 
outflow, and M is the Mach number of the flow relative to 
the ambient material. 

The best-flt from flgure 2 suggests tan(0) = a ~ 0.25. 
Typical values for a momentum dominated outflow might 
be M ~ 10, and p/p outflow 

~ 10^ which would give 
tan{9) ~ 0.25, using equation (|3ip . In the buoyancy domi- 
nated regime, the flow speed can be described by M ~ 1, 
for which equation (|3ip suggests p/poutflow ~ 10. The num- 
bers used in this calculation are likely to be common to 
all clusters, suggesting that entrainment of material might 
well be what sets the effective opening angle. Furthermore, 
if the outflows can be described as poutflowW^ = constant, 
as expected for momentum conservation, then the effective 
opening angle will also be constant, according to equation 

In flgure 2, A478 appears not to be well-described by the 
best-fit through the data. This may be because the best-fit 
assumes a uniform effective opening angle for all the outflows 
in the cluster sample, and the outflows in A478 typically 
have a larger effective opening angle. There does not seem 
to be any particular reason to expect this ; altough A478 is 
the most massive cluster in this sample, it does not exhibit 



a gravitational acceleration that is significantly different in 
magnitude to any other clusters in the sample. 

As an alternative explanation, it is interesting to note 
the la rge temperature decrement in the cluster centre of 
A478 (|Sun et al.l 120031 ') , which is often taken to be a sig- 
nature of a strong inflow. The AGN also appears to be in- 
jecting energy at a very low rate compared to the radiative 
losses (e.g. see iBirzan et al.ll2004l ). As a result, it might be 
no coincidence that A478 seems to have a larger mass within 
^■€001 than would be expected from the best-flt to the data. 

Table 1 lists the X-ray luminosities, weights and veloc- 
ity dispersions used to plot flgure 2. For consistency with 
our earlier calculations it was assumed that = k\,T / prUp, 
where the temperature was taken at its highest value within 
the cooling radius. The values in table 1 have been rounded 
to two signiflcant flgures. For this reason, and the small sam- 
ple size, the uncertainties were not calculated since they do 
not provide any extra information. 

An alternative explanation for this correlation is not 
obvious. So, it seems plausible to conclude that energy in- 
jection by AGN does reshape the gas profiles in the central 
regions of clusters. As a consequence, it is not a surprise that 
heating and cooling rates match relatively well - because the 
systems were reconfigured into such a state by the AGN. 

Further insight into the interaction between the AGN 
and the ambient material can be gained by studying equa- 
tion (|29[) in more detail. If the equilibrium condition is met 
at all radii, then 

Lx(< R) = - {^^ (^^) aM{R)g{R). (32) 

The X-ray luminosity is defined by Lx ~ / n^A(r)dl/, 
where n is the electron number density of the gas, and A(r) 
is the temperature dependent cooling function. For typical 
cluster temperatures we expect the emission to be domi- 
nated by thermal bremsstrahlung so A(T) ~ T^^'^ . 

Differentiating both sides of equation (I32p . with repsect 
to volume, leaves 



n^A(r) = - 

Factorising g{R)p{R), we find 
n^A(r) = 



g{R)p{R) + 



Mt{R) dg{R) 



47ri?2 dR 



(33) 



'-l)(§;^]agiR)piR) 



1 + 



din g{R) 



dlnMT(i?) 



(34) 

The gas distribution in clusters is usually approximated by 
a /3-profile, which is written 

po 



p(i?) = 



[1 + [R/ro] 



21/3 ■ 



(35) 



where po is the central density and ro is the scale height. 
Assuming the the gas temperature is independent of radius, 
then Vp — {khT / pmp)dp/dR, so the gravitational accerla- 
tion is 



Vp(-R) 

p(R) 



kbT 2/37? 



prup 



i?2)' 



(36) 



and equation (|34p can be rearranged to obtain the critical 
cooling time 



cooljCrit — 



5 nfcbT 



2 n^A{T) 



1 

3/3^ 



r-g + R' 
R 



(37) 
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Table 1. Table showing the data used to plot figure 2. Column 1 lists the names of the clusters, column 
2 contains the observed X-ray luminosities within the cooling radius, Lx(< J'cooOi column 3 gives the 
(''cool )9 ('"cool) values, derived from temperature and density fits to the data. Column 4 lists the velocity 
dispersions defined as = k^^T/ finip. 



Name 


ix(<rcool)/104'ergs-i 


-MT(»'cool)5(»'cool)/10^^gcms 2 


cr/kms ^ 


Virgo 


9.8 


0.55 


620 


Perseus 


670 


20 


1100 


Hydra 


250 


11 


780 


A2597 


430 


18 


800 


A2199 


150 


3.7 


830 


A1795 


490 


23 


970 


A478 


1220 


62 


1100 



8x10^' 



6x10™ - 



4x10 - 



2x10™- 




2.0x10™ 4.0x10™ 6.0x10™ 8.0x10™ 1.0x10-" 1.2x10'' 
observed X — ray luminosity/ velocity dispersion (g cm s"^) 



Figure 2. Plot for equation ((29J using data for the Virgo, Perseus, Hydra, A2199, A2597, A2199, A1795, A478 taken from Pope, et al. 
(2006) and references therein. The plot confirms that Mt(?'cooi)9(''cooi) is linearly related to the X-ray luminosity within the cooling, as 
the model predicts. The solid line shows the best-fit to the data, which yields an opening angle for the outflow of f2/47r ~ 0.028. 



The minimum occurs at R - 
simplifies to 



ro, at which point the function 



tcool.cHt ~ ( ^) (5 



(38) 



This relation is only approximate since it was assumed that 
the term in square brackets in equation (|34|l is 1. How- 
ever, the approximation is reasonable because we expect 
d In g/d In Mt ^ near ro . 

For cooling times below this threshold, material can 
build up in the core of a cluster, because the energy in- 
jection rate is insufficient to drive it away. This build up of 
material is associated with the star formation observed in 
clusters with short central cooling times. 

As a general example, let ro — 30kpc, a = 10^ kms~^, 
P = 3/4 and r2/47r = 0.028, so that ^cooi.c rit ~ 0-5 Gy r This 
agrees well with the value found by Ra ffertv et al.l ()2008l ) 
below which star formation appears to 'switch on'. Note 
that this is also similar to the entropy threshold found by 
ICavaenolo et all (|2008l i. 

In table 2, we list the critical cooling times calculated 
for the Virgo, Perseus and Hydra clusters. These objects 
were chosen because /3-profiles are used to fit the gas den- 
sity profiles, and the clusters are sufficiently resolved that 
the fits are reliable to small radii, ro and /3 are taken from 



Pope et al.l ||200^. and references therein (see [P avid et al.l 
200 ll : I Ghizzardi et al ] |2004l : ISanders et al.l[20o3 ). The Virgo 
and Perse us clusters ca n bo describe d using a doub le j3- 
profile fsee lChurazov et a l. 2003; Ghiz zardi et al.|[2004h with 
different scale heights and /3s. For the Virgo cluster, we took 
the parameters for the profile with the larger scale height, 
and used a single /3-profile fit for the Perseus cluster. 

Simple analysis also suggests that the ratio ro /cr should 
be fairly constant. According to the virial relations we expect 
oc M/ro oc ro, so that ro/a should be approximately 
constant. As a result, a cooling time of 0.5 Gyr should mark 
a transition in behaviour in all clusters, as appears to be the 
case. Whenever, the cooling time of the gas, drops below this 
critical value, we should expect to observe the consequences 
of material building up in the centre of a cluster. 



It has b een suggested by ICavagnolo et al.l (|2008l i and 

ISokeJ i 20081 ) that the presence of the critical cooling time 
is a possible indication th at central AGN are fuell ed by the 
cold feedback mechanism ijPizzolato fc Sokeijl2005h . Instead, 
the model presented here suggests that the phenomenon oc- 
curs directly because of the buoyancy injected by the AGN. 
This places no requirements on the AGN fuelling mecha- 
nism, except that it must balance the radiative losses within 
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Table 2. Table showing the data for the critical cooling time. Column 1 gives the cluster name, column 2 
lists the core radius of the density profile, column 3 gives the /3 value of the density profile and column 4 
shows the critical cooling time according to equation (1371 . using the appropriate a values from table 1. The 
opening angle was taken to be Q/4n = 0.028, from the best-fit to equation l|29|l . The values of ro and /3 are 
taken from Pope, et al.(2006), and references therein. 



Name 


ro/kpc 


/9 


tcooLcrit/10^ yrs 


Virgo 


23.3 


0.71 


0.5 


Perseus 


28.5 


0.81 


0.3 


Hydra 


18.6 


0.72 


0.4 



the cooling radius. Cold feedback still provides a possible 
mechanism for achieving this. 

An alternative view, presented bv lVoit et al.l (|2008l ), is 
that the onset of star formation and optical emission oc- 
cur below a critical entropy. They showed that for a cen- 
tral entropy greater than Ko ~ SOkeVcm'^, thermal conduc- 
tion transports energy at a rate greater than can be radi- 
ated away by the hot gas, assuming that the conduction is 
not heavily suppress ed. In this case entropy is defined as 
K = ki^Tn-^^^ (e.g. iLlovd -Davies et al where fcb is 

the Boltzmann constant, T is the temperature and rie is the 
electron number density. Material below the critical entropy 
would be thermally unstable, able to cool, and could poten- 
tially fuel the central AGN. This expla n ation fits neatly into 
the framework outlined bv lGuo et al.l (|2008l ) in which non 
cool-core clusters can be stabilised by thermal conduction, 
whereeis cool core clusters require AGN feedback and ther- 
mal conduction for stability. Their view is not inconsistent 
with the model presented here. For example, the model dis- 
cussed here focusses on clusters in which AGN feedback is 
important, and thermal conduction cannot balance radia- 
tive cooling. However, complete consistency would require 
that AGN heating either i) does not raise the central en- 
tropy above ~ 30keVcm^, otherwise the gas would never 
be able to cool to fuel the AGN again, or ii) that thermal 
conduction is heavily suppressed. 

One difficu lty with the possibility suggested by 
IVoit et al.l ()2008l ) is that the value of the thermal conduc- 
tion suppression factor is highly uncertain. In contrast, the 
model presented here predicts a critical cooling time which 
contains no free parameters. Therefore, probably the best 
way to distinguish between the models would be to deter- 
mine whether, and by how much, thermal conduction is sup- 
pressed. Unfortunately, this is a question that cannot easily 
be answered. 

If the model presented here does describe some of the 
real mechanisms at work in clusters, perhaps the most im- 
portant aspect to highlight is the self-tuning ability of AGN 
feedback. This is true in both elliptical galaxies and clusters. 
The tuning is achieved by slightly different mechanisms in 
each case. In elliptical galaxies, an equilibrium is reached 
when the AGN has ejected so much gas that its fuel sup- 
ply becomes limited and it can no longer exceed the critical 
threshold required to eject gas from the gravitational poten- 
tial. In clusters, the AGN injects energy at a particular rate 
given by the details of the accretion process. This energy in- 
jection has the twin effect of heating the gas and also moving 
it about within the cluster potential. Within the radius that 
the AGN heating balances the radiative losses, the heating 



also controls the weight of cooling material, and therefore 
its own fuel supply. 



5 SUMMARY 

We have constructed an analytical model to investigate the 
observable effects of kinetic and thermal AGN feedback on 
elliptical galaxies and galaxy clusters. The results show that 
there is a critical momentum injection rate above which an 
AGN outflow will eject mass from an elliptical galaxy. For 
an AGN fuelled at a rate proportional to the classical mass 
cooling rate, this leads to Lx oc o-^° which agrees well with 
observations fe.g. iMahdavi fc Geller 20pi|) and reasonabl y 
with simulations (e.g. iBower et al.ll2008l : lOave et al.ll2008l ). 
If the AGN fuelling rate is only proportional to Lx this 
would lead to Lx oc cr*. Presumeably these represent the 
extremes, so that a different fuelling mechanism might lead 
to an intermediate scaling. 

In more massive environments, the work done by ther- 
mal energy injected by the AGN also becomes significant. 
This mechanism should lead to Lx oc cr^, assuming the AGN 
fuelling rate is proportional to the classical mass cooling 
rate. It is not clear whether the observations confirm this, 
due to the large scatter. Again, if the AGN fuelling rate is 
oc Lx, this would lead to Lx oc a^. Scalings between and 
a"^ are also, no doubt, possible in this regime, depending on 
the details of the AGN fuelling mechanism. 

Above a ~ 500kms~^ AGN heating of either sort is 
unable to affect the global Lx — f relation. This corresponds 
to a temperature of roughly 1-2 keV an d can explain th e 
observed break in the Lx — L relation (e.g. lXue fc Wtj|2000l ). 
In order for this break to be noticeable, it is required that 
the Lx — cr relation is steeper for ellipticals than for clu sters. 
Since Lx oc a* '^ in clusters (iMahdavi fc GellerllioOll ). this 
implies that the AGN fuelling rate probably does scales more 
like Lx/f^ than Lx. 

In clusters of galaxies, the stronger gravity means that 
outflows become buoyancy dominated above a certain length 
scale. Simple estimates suggest that ~ 20 kpc might be typ- 
ical. It has also be shown that buoyancy is equivalent to 
the work done by thermal energy injection and since typical 
cooling radii in clusters are ~ 100 kpc, global considerations 
of AGN heating in clusters can be treated as if the energy 
injection is purely thermal. 

It is not possible for an AGN to completely eject mate- 
rial from a cluster, however, it is possible for the energy in- 
jection to redistribute material within the central regions of 
the cluster. This can lead to the self-tuning of AGN feedback 
so that heating roughly matches cooling within the cooling 
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radius. A direct consequence of this is that the mass of hot 
material within the cooling radius should be proportional to 
Lx{< rcooi)/[o"t;(rcooi)]- This appears to be true, at least for 
a sample of seven nearby clusters. The only possible excep- 
tion is A478, which also displays evidence of a strong inflow 
and minimal AGN heating and might therefore be expected 
to contain excess material in the central regions. Further- 
more, the same physical processes lead to the derivation of 
a critical cooling time of ~ 0.5 Gyr, below which material 
can pile up in the cluster centre. This is very close to the 
observationall y determined u n iversa l cooling time of ~ 0.5 
Gyr, found by iRaffertv et af] (|2008l ). below which star for- 
mation and AGN activity appear to be triggered. 
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